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Abstract 

A computational study of the wing for the distributed electric propulsion X-57 Maxwell airplane configuration at cruise 
and takeoff/landing conditions was completed. Two unstructured-mesh, Navier-Stokes computational fluid dynamics 
methods, FUN3D and USM3D, were used to predict the wing performance. The goal of the X-57 wing and distributed 
electric propulsion system design was to meet or exceed the required lift coefficient 3.95 for a stall speed of 58 knots, 
with a cruise speed of 150 knots at an altitude of 8,000 ft. The X-57 Maxwell airplane was designed with a small, high 
aspect ratio cruise wing that was designed for a high cruise lift coefficient (0.75) at angle of attack of 0°. The cruise 
propulsors at the wingtip rotate counter to the wingtip vortex and reduce induced drag by 7.5 percent at an angle of 
attack of 0.6°. The unblown maximum lift coefficient of the high-lift wing (with the 30° flap setting) is 2.439. The stall 
speed goal performance metric was confirmed with a blown wing computed effective lift coefficient of 4.202. The lift 
augmentation from the high-lift, distributed electric propulsion system is 1.7. The predicted cruise wing drag coefficient 
of 0.02191 is 0.00076 above the drag allotted for the wing in the original estimate. However, the predicted drag overage 
for the wing would only use 10.1 percent of the original estimated drag margin, which is 0.00749. 


Nomenclature 

Cp drag coefficient Viratio ratio of tip speed to freestream velocity 
Cp HLN drag coefficient, high-lift nacelles contribution = W aircraft weight, lb 
Cp pylons drag coefficient, pylons contribution y axis along the wing span, in. 
Cp.rw drag coefficient, wingtip nacelles contribution = y* nondimensional first node height in boundary layer 
Cp,wing drag coefficient, wing contribution Veo nondimensional first cell centroid height in boundary layer 
Ce skin friction coefficient 
C, lift coefficient Symbols 
cl sectional lift coefficient a angle of attack, degrees 
Cerf effective lift coefficient: C,+ Cy prop A delta 
Cromae maximum lift coefficient p density 
Ch prop lift coefficient from the contribution of Acronyms 

propeller thrust in lift direction BSL Menter k-w basic turbulence model 
Cm pitching moment coefficient CFL pseudo time advancement Courant-Friedrichs-Lewy 
Cy pressure coefficient DEP distributed electric propulsion 
Che f reference chord, in. HLN high-lift nacelles, including pylons 
Cr thrust coefficient HP horse power 
Co torque coefficient KCAS knots calibrated airspeed 
D drag force KEAS knots equivalent airspeed 
d propeller diameter, ft. KTAS knots true airspeed 
h altitude, ft. LM Langtry-Menter transition model 
Kr normalized thrust coefficient mph miles per hour 
Kg normalized torque coefficient QCR quadratic constitutive relation 
M Mach number RPM revolutions per minute 
P pressure, Ibf/in? SA Spalart-Almaras one equation turbulence model 
q dynamic pressure SARC SA rotation and curvature correction 
Re Reynolds number based on C,¢¢ SCEPTOR Scalable Convergent Electric Propulsion Technology and 
Ss wing reference area, ft Operations Research 
T temperature, °F SST Menter’s Shear Stress Transport model 
Vv freestream velocity, ft/sec TN wingtip nacelles 
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I. Introduction 

The NASA New Aviation Horizons Initiative has a goal of green aviation to alleviate the negative 
impacts associated with the expected doubling in growth of passenger trips per year from now to 2035, and 
to take advantage of the economic potential of being aeronautical leaders in the future aviation growth [1]. 
Clean, quiet, and even supersonic airplanes are planned within the NASA X-plane research approach to 
improve energy efficiency, and develop technologies and new propulsion systems that will enable 
innovative airplane designs. The X-57 Maxwell, or Scalable Convergent Electric Propulsion Technology 
and Operations Research (SCEPTOR) airplane is one of several airplanes worthy of the X-plane status. One 
goal of the all-electric X-57 Maxwell technology demonstrator is to prove that significantly reducing wing 
area (to improve cruise efficiency) can be done without compromising takeoff and landing performance. 
This is accomplished by application of distributed electric propulsion (DEP) to effectively blow a higher “q” 
(dynamic pressure) over the wing and flap during takeoff and landing conditions. 

The new DEP wing system for the X-57 airplane demonstrator will be installed on a Tecnam P2006T 
aircraft (fig. 1) by removing the original wing and installing the cruise-optimized DEP wing (fig. 2). The 
main performance goals are a 58 knots equivalent airspeed (KEAS) stall speed (scaled from the original 
P2006T gross weight and stall performance), at 3,000 Ib gross weight equating to a C, = 3.95, while 
achieving a cruise speed of 150 knots true airspeed (KTAS) at 8,000 feet. It is projected that the all-electric 
X-57 Maxwell airplane will achieve 5 times lower energy use than the original P2006T. The efficiency gains 
will be accomplished with aerodynamic improvements and by the increased propulsion system efficiency of 
the electric motor/battery system. The sizing design study of the wing presented in Reference 2 resulted in a 
wing design with a wing loading of 45 Ib/ft?, a wing area of 66.67 ft’, an aspect ratio of 15, and a cruise C,= 
0.75. The higher aspect ratio for the new wing is needed to minimize induced drag at cruise, since the new 
wing’s cruise lift coefficient is much higher than the original P2006T. The original P2006T has a wing 
loading of 16.365 lb/ft?, a wing area of 158.88 ft”, an aspect ratio of 8.8, and a cruise C, = 0.275. The 
specially designed X-57 airfoil is tailored for a cruise lift coefficient of 0.90 and incorporates a 25% chord 
flap. The flap design uses a single-pivot displaced hinge with a 30° maximum deflection. 


Il. Configurations 

The configurations analyzed in this computational study of the wing for the DEP X-57 Maxwell 
airplane were the cruise wing and the high-lift wing, which has a 30° flap setting (figure 3). The X-57 
Maxwell wing has been through several design improvement revisions and the configurations that are 
presented in this paper are the 4.1 design revision. The X-57 Maxwell airplane has 12 high-lift propellers 
mounted on nacelles upstream of the wing leading edge that are positioned in an alternating fore- and aft- 
staggered pattern. The high-lift propellers are designed to fold smoothly onto the nacelle for the cruise 
configuration (fig. 3(a)), which is referred to as the cruise wing in this paper. The high-lift propellers (not 
shown) are positioned on the nacelles in an alternating fore- and aft-staggered pattern for the high-lift blown 
wing (30° flap shown in red in fig. 3(b)). The cruise wing without HLN (fig. 3(c)) and the high-lift wing 
without HLN (fig. 3(d)), were also analyzed to investigate the impact of the pylon and nacelle geometry on 
wing performance (lift coefficient discussed in this paper). The fuselage and tail were not included in this 
study but will be included in future computational fluid dynamics (CFD) methods analysis of the X-57 
airplane. 


III. High-Fidelity Distributed Propulsion Analysis Tools 

Two unstructured-mesh, Navier-Stokes CFD methods were used to define the wing performance, 
FUN3D [3] and USM3D [4, 5]. Additionally, wing performance was compared with other computational 
methods (OVERFLOW and STARCCM+) in reference 6, with the goal to characterize and provide further 
confidence in the computational results. 

FUN3D is a node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. 
Fully turbulent predictions were computed with FUN3D using the Spalart-Allmaras one-equation turbulence 
model with rotation and curvature correction (SARC) [7] and the quadratic constitutive relation (QCR) [8] of 
Spalart. Additionally, some cases were also computed with two-equation turbulence models; Menter k-w 
basic (BSL) and Shear-Stress Transport (SST) [9]. For this work, FUN3D was the primary method used to 
define Cmax for the cruise wing and the high-lift wing with a 30° flap. All of the FUN3D simulations for 
this study used an actuator disk representation for the high-lift propellers, with thrust and torque coefficient 
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data being inputs to the solution. The thrust and torque coefficient data were derived from the XROTOR 
[10] blade element momentum analyses. 

USMS3D is a cell-centered, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. 
Three turbulence models were used in this study; SA [7-8], SST [9], and the Langtry-Menter (LM) transition 
model [11-12]. In order to understand the drag benefit that might be possible with laminar flow, the USM3D 
flow solver with the LM transition model was used for predicting transitional flow along the wing. Fully 
turbulent solutions were also computed with the SA and SST turbulence models, using both USM3D and 
FUN3D, to provide a baseline for comparison to the laminar flow prediction. In this paper, USM3D was not 
used for powered solutions. 


A. Grid Generation 

GridTool/VGRID [13] software was used to generate the meshes for FUN3D and USM3D simulations. 
GridTool is used to take a supplied geometry file and create the necessary points, curves, patches and grid 
topology (sources) to define the surface for grid generation. VGRID is a tetrahedral unstructured grid 
program that can either be run interactively or in a batch mode with grid spacing corresponding to the 
strength of user-defined sources placed in the domain. The advancing layers method [14] is used to generate 
the boundary layer and the advancing front is used to generate the farfield tetrahedral mesh. POSTGRID is 
used for additional post-processing, to close any open pockets that do not have cells and to improve the grid 
quality. 

An unstructured, mixed element mesh was used for the FUN3D simulations. The mixed element mesh 
had prisms in the boundary layer and tetrahedral cells in the farfield, with pyramids to transition between the 
boundary layer and farfield cells. Semispan wings were used in the computational model since all 
calculations were computed at 0° angle of sideslip and little lateral variation was expected. A full span 
solution was computed for the blown, high-lift wing (30° flap) with high-lift nacelles at a = 0°. The 
differences between the full span and semispan configurations for lift, drag, and pitching moment coefficient 
were 0.2%, 0.8% and 0.1%, respectively. These small differences between the solutions confirmed the 
approach of using a semispan mesh. 

For all grids, a flat-plate, turbulent boundary layer calculator was used with the reference chord length, a 
Reynolds number based on reference chord length, the intended y*, and an intended number of layers of 
cells within the boundary layer, to determine the first node height and expansion rate of the grid for input to 
VGRID. Previous results from the LEAPTech wing (not shown in this paper) showed that using an intended 
y*= 0.67 for a medium mesh had solution y* values up to 2 for the blown, high-lift wing solutions, even 
though the solution yt were less than | for the unblown wing solutions. Values of solution y* that are less 
than | are recommended to adequately model the boundary layer flow. 

For all grids, the farfield boundaries were extended to approximately 100 chord lengths away from the 
wing and the reference chord length was C;e¢= 25.56 in. 

For the FUN3D code, the semispan mesh of the cruise wing with wingtip nacelles had 43.9 million 
mesh points and the semispan mesh for the cruise wing with high-lift and wingtip nacelles had 47.9 million 
mesh points. As recommended in reference 15, both meshes were generated with an intended y* = 0.3 for an 
extra fine mesh, and 32 grid layers were specified within the boundary layer. The corresponding VGRID 
spacing for the first node height was 6.3e-5 in. as determined by the boundary layer calculator for the cruise 
mesh at Re = 2.83 million, and 150 KTAS. 

The same grid created for the cruise wing with high-lift and wingtip nacelles was used for USM3D. The 
number of cells are important for USM3D since it is a cell-centered code, and the 47.9 million mesh points 
grid had 282.5 million cells. This means that for the same grid USM3D has 6 times more resolution than 
FUN3D. Consequently, a smaller grid could be used for USM3D, however, since time was limited, the same 
grid was used and another grid with similar resolution as the FUN3D solution was not made. The advancing 
layers method forms prismatic layers that are divided into three tetrahedral cells. The distance from the 
surface to the centroids of the tetrahedral boundary cells is one-fourth that of the first layer of nodes, 
resulting in the yé, = 0.25 - y* [16]. Therefore, the FUN3D mesh used for running the USM3D code has 
the yé, = 0.075. 

The FUN3D semispan mesh for the high-lift wing (30° flap setting) and tip nacelle had 149 million 
mesh points. Adding the high-lift nacelles brought the grid to 153 million mesh points. Both meshes were 
generated with an intended y* = 0.3 and 32 grid layers within the boundary layer. The corresponding 
VGRID spacing for the first node height was 1.3e-4 in. for the high-lift wing mesh at Re = 1.33 million, and 
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58 KTAS. 

No special patching is required within the grid for FUN3D to simulate the propellers, except clustering 
of cells to resolve the flow. The clustering of cells was achieved using VGRID volume sources during grid 
generation. Grids for powered runs with the USM3D code do require patches to create the actuator disk, 
which simulate the propellers, however, powered solutions were not computed with USM3D for this paper. 


B. Computational scheme 

The process for running the FUN3D code for all of the cases shown was to start the solutions as steady- 
state with a local time stepping method. If needed, solutions were switched to time-accurate, global time 
stepping, to obtain acceptable convergence. In nearly all of the cases with separated flow, which were near 
the stall condition, the final converged solution was computed with the time-accurate global time stepping 
method. 

Since the USM3D code was run on the cruise wing at cruise conditions where little flow separation was 
expected, the solutions were obtained using the steady-state, local time stepping method and converged well. 


C. Boundary Conditions and Actuator Disk Modeling 

A symmetry plane boundary condition is used for the semispan grids while the farfield boundaries 
utilized the Riemann invariants boundary condition. All of the solid surfaces were set with the no-slip 
boundary condition. 

An actuator boundary condition was used to simulate the high-lift (HL) propellers. There are 6 options 
to loading the rotor in FUN3D but the one implemented for this work was a body force based on the optimal 
distribution of Goldstein [17] implemented as described by Stern, Kim, and Patel [18]. Inputs for the 
actuator boundary condition include tip radius (11.34 in.), hub radius (3.06 in.), and the (x, y, z) location of 
the hub centers for each high-lift propeller. Additionally, the Vt,g;i,, the normalized thrust coefficient and 
the normalized torque coefficient were required for each propulsor. The Vt,q;;,, as defined by equation 1 
was calculated to be equal to 4.5976. Positive values of Vt,g¢ig results in the propellers on the right wing to 
rotate in a clockwise direction from the pilot’s view, a rotation that opposes the wingtip vortex. The torque 
and thrust coefficients, equations 2 and 3, respectively, were calculated using 4548 RPM, a torque of 15.9 ft- 
Ibf, a thrust of 49.7 Ibf, and the propeller diameter (d) of 1.89 ft. The normalized torque coefficient (Kg = 
0.0125) and normalized thrust coefficient (Kr = 0.0368) were computed with equations 4 and 5 and used in 
the actuator boundary condition input file. 


Vtratio = [x Ps (=) Ps d| /V (1) 
Co = Torque/[p(-)"4°] (2) 
Cy = Thrust /[pCj)°*] (3) 
Ky = = Gs (4) 
Ky = Cr (5) 


D. Flow Conditions 
Computational solutions of the cruise wing were computed at a cruise speed of 150 KTAS and at an 
altitude of 8,000 ft. Computational solutions of the high-lift wing (with a 30° flap setting) were computed for 
a landing speed of 58 KTAS at sea level. 


E. Convergence 
The criteria used to monitor and determine solution convergence was a drop of at least two orders of 
magnitude for the flow solution residual. In addition, for an interval of 1,000 iterations, it was required that 
for steady-state solutions the standard deviation of C;, be less than 0.0012 and Cp be less than 0.0020. 
The FUN3D steady-state solutions for the cruise wing with tip nacelle all converged to a standard 
deviation less than 0.0001 for both C,, Cp, and C,,. Steady-state cases for the cruise wing with high-lift and 
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wingtip nacelles (a<10°) converged to a standard deviation less than 0.0005, 0.00006, and 0.0002 for C,, 
Cp, and C,,, respectively. 

The USM3D code was run on the cruise wing at an angle of attack of 0.6°. The USM3D SA solution 
converged with standard deviations of 0.00009, 0.000003, and 0.00002 in C,, Cp, and C,,, respectively, over 
2000 iterations. The USM3D SST solution converged with standard deviations of 0.00007, 0.0000004, and 
0.00003 in C,, Cp, and C,,, respectively, over 2000 iterations. The USM3D LM solution converged with 
standard deviation of 0.00035, 0.00003, and 0.00014 in C,, Cp, and C,,, respectively, over 6000 iterations. 


IV. Results 

The traditional lift coefficient (C,) and the effective lift coefficient (C,.¢7) are discussed throughout the 
paper. The effective lift coefficient includes the addition of the vectors of propeller thrust and normal force 
in the lift direction, to the lift coefficient computed from the CFD pressures. The propeller contributions 
have a small effect at low angles of attack, but can be on the order of AC, = 0.10 to 0.20 at a = 10° in the 
high-lift DEP cases. 

A component drag breakdown was conducted in the NASA SCEPTOR DEP demonstrator (X-57 
Maxwell) sizing study [1] and the pie chart of the drag elements is presented in Figure 4. The drag 
components are summed up in Table | and the resulting airplane drag coefficient required to meet the cruise 
speed goal is Cp = 0.05423. Note Cp is nondimensionalized with the new X-57 wing area, S = 66.67 ft?. A 
drag margin of Cp = 0.00749 was built into the airplane drag coefficient required to meet the cruise speed 
goal, to account for uncertainties in the estimate. The drag margin is about 13.8% of the total estimated 
airplane drag. The drag components attributed to the wing are summed up in Table 2. The drag of the wing 
(induced + friction + profile + cruise nacelles + high-lift nacelles) should not exceed Cp = 0.02115. 

The FUN3D results of C;, and Cp for the cruise wing with the high-lift nacelles are shown in Figure 5. 
The unpowered cruise wing achieves C; ma, = 1.7 at a = 16° (fig. 5(a)). A high lift coefficient is needed at 
the cruise condition because the DEP wing is designed to be small with a high aspect ratio for efficiency. At 
an angle of attack of a = 0.6°, the unpowered wing has a lift coefficient C, = 0.756. When cruise power is 
included in the solution, the lift coefficient increases to C, = 0.780. The cruise propulsors at the wingtip 
rotate counter to the wingtip vortex and do provide some induced drag reduction (fig. 5(b)). The cruise 
propulsors reduce drag by 16 counts at a = -2°, by 23 counts at 0° and 2° angles of attack, and by 24 counts 
at 0.6° angle of attack. At a = 0.6°, the 24-count drag reduction from the wingtip propulsor accounts for a 
total wing drag reduction of 7.5 percent. The drag coefficients for these fully turbulent flow solutions at a = 
0.6° is Cp = 0.03254 for the unpowered case and Cp = 0.0301 for the cruise power case because of the 
induced drag reduction. The design target wing drag coefficient needed for the cruise condition is Cp = 
0.02115. The CFD turbulent drag coefficient does include the drag on the wing planform that is inside the 
fuselage. It also does not account for achievement of laminar flow on the wing. These corrections are 
explained as follows. The CFD semispan grid extends from the fuselage centerline (y = 0.0 in.) to the 
wingtip, however, the half-width of the fuselage extends to the wing span location of y = 24.0 in. The 
portion of the new DEP wing that contributes to aircraft drag will be outside of station y = 24 in. The 
spanwise section drag coefficient (cz), from the cruise wing CFD solution, was first integrated from the 
centerline to the wingtip. Then the spanwise section drag coefficient (ca) was integrated from wing station y 
= 24 in. to the wingtip. The difference in the two totals of the drag coefficient, A Cp = -0.00480, is the drag 
computed on the wing planform that is inside the fuselage. This yields a fully turbulent, wing drag 
coefficient that contributes to aircraft drag of Cpy>24= 0.03010 — 0.00480 = 0.02530. 

However, the drag breakdown of the performance sizing assessment assumed that some laminar flow 
would be achieved on the wing. An analysis was initiated to determine the drag savings due to laminar flow 
on the wing that might be expected in flight on the X-57 airplane, compared to the fully turbulent predictions 
from FUN3D. The coauthors were unable to get the FUN3D turbulence model with transition to work, using 
several grids made with VGRID, POINTWISE and GridEX. In order to understand the drag benefit of 
laminar flow that may be present at cruise, a grid of the cruise wing configuration with high-lift nacelles was 
run with the USM3D flow solver with several turbulence models and the LM transition model. Two 
turbulence models were run with USM3D as a comparison for the transition model approach and as a check 
against the existing FUN3D results. The analysis at a = 0.6° was for the unpowered cruise airspeed of 150 
KTAS, T = 30.5°F, h = 8000 ft, and Re = 2.83 million. Figure 6 shows the skin friction coefficient (Cy) on 
the upper and lower surfaces of the wing for the cases analyzed. Figure 6(a) shows the extent of the laminar 
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flow predicted with USM3D LM, with low values of C; to approximately 70 percent chord on the upper 
surface and 60 percent chord on the lower surface, where Cy rises to values around 0.005 on the upper 
surface and 0.003 on the lower surface, as the flow transitions to fully turbulent flow. In contrast, the fully 
turbulent flow cases (figs. 6(b) to 6(e)) have high Cy values near 0.005 to approximately 70 percent chord 
and the contours are relatively similar between the solutions on the upper surface. However, there are more 
differences in Cy on the lower surface between the fully turbulent cases, with similar contours for the 
FUN3D solutions, while USM3D predicts higher Cy values for SA than FUN3D, and even higher values of 
C; for SST than SA. At this angle of attack a@ = 0.6°, the lower surface flow is more influenced by the flow 
around the pylons and nacelles than the upper surface. Table 3 shows the total drag coefficients and the 
contributions of the wing, wingtip nacelles, high-lift nacelles and pylons to the total drag for the fully 
turbulent and transitional flow solutions from FUN3D and USM3D. The total drag coefficient for the fully 
turbulent models all vary between Cp = 0.03254 and Cp = 0.03557. The drag coefficient of USM3D with 
the LM transition turbulence model is lower, Cp = 0.02938, with laminar flow on the wing. The LM 
transition turbulence model does predict laminar flow on the high-lift nacelles as well. However, because of 
the intended folding propeller installation, practical laminar flow is not expected on the HLN. Therefore, 
when only the wing drag is considered, and the same code (USM3D) is compared to reduce uncertainty, then 
the resulting drag reduction due to laminar flow achieved on the base wing planform is a drag reduction of 
ACp = -0.00397. This drag reduction with laminar flow is compared to the fully turbulent result from the 
SST turbulence model as the basis for comparison. The LM transition model is highly validated against the 
SST turbulence model [12, 19-20]. 

Note the computed benefit of laminar flow, A Cp = -0.00397, includes the wing from the centerline to 
the wingtip (CFD grid). However, as stated before, the wing planform inside of fuselage station y = 24” is 
within the fuselage and does not contribute to the aircraft drag. The ratio of the wing planform area from the 
wing centerline to station y = 24 in., to the total wing area is 0.146. Therefore, if we only count the drag 
benefit of laminar flow on the wing outside of the fuselage, then laminar flow drag reduction is calculated as 
ACD laminar.y>24 = -0.00397 * (1 - 0.146) = -0.00339. The final computed wing drag coefficient is Co wing = 
Cp,y>24 + ACD laminary>24 = 9.02530 — 0.00339 = 0.02191. The target estimate for the wing drag is 0.02115 + 
0.00749 = 0.02864, which includes 74.9 counts of drag as margin. Thus, only 10% of the originally planned 
drag margin is required to meet the design cruise goal. 

CFD analysis of the high-lift wing are examined next. The traditional and effective lift coefficients for 
the blown and unblown high-lift wing at a freestream velocity of 58 KTAS are presented in figure 7. This 
figure illustrates the increase in lift coefficient over the range of angle of attack from the unblown wing 
(purple curve), to the blown wing (red curve), as well as the extra benefit to lift coefficient from the thrust 
vector (blue curve). The unblown C; ma, of the high-lift wing is 2.439. The blown high-lift wing, with the 
DEP system operating at 164.4 hp (13.7 hp/prop), achieves an effective C, max = 4.202 at a = 12°. The 
traditional, maximum lift coefficient is 3.89 at a@ = 12°, without the addition of the propeller vectors 
contributing to lift. At a = 10°, the blown wing C, is 4.15 and the unblown wing C, is 2.4, a lift 
augmentation of 1.7. The impact of the lift augmentation is illustrated in figure 8, a comparison of the 
spanloading between the unblown and blown high-lift wing at 58 KTAS and @ = 10°. The pink vertical lines 
are placed at the propeller tips to represent the propeller diameter, and the green vertical lines are positioned 
at the center of the prop hub. The blue curve (s4.1g4m088a10p0.sect_forces) is the spanwise sectional lift 
coefficient (c;) distribution for the unblown high-lift wing with HLN, and the slight jaggedness in the curve 
at the prop hubs represents the effect adjusting the chord that included the nacelles to a wing chord along the 
span. The red curve (s4.1g4m088a1l0p2.sect_forces) for the c; of the blown high-lift wing with HLN is 
substantially higher than the unblown high-lift wing with the peaks in c; at the prop upwash side of the 
nacelle. Local spanwise section lift coefficients as high as 4.5 to 5.4 are achieved with the high-lift blowing. 
It can also be seen that the blowing increases the local section c; outside of the blown region. The lift 
augmentation due to high-lift blowing is represented by the area between the blue and red curves. 

Figure 9 shows the difference between two turbulence models on the computed effective lift coefficient 
for the blown, high-lift wing. In general, the turbulence models match fairly well through a = 8°, and the 
BSL model predicts a later stall angle of attack by at least 3°. The BSL two-equation turbulence model also 
predicts a higher effective lift coefficient (C, ma, = 4.444), at a higher angle of attack (a = 14°), than the 
SARC+QCR turbulence model previously discussed. Although the CFD solutions have different 
C,,max Predictions, the important take away is that even the conservative estimate is exceeding the design 
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C,max and predicts that the stall speed goal can be met. The advanced one-equation SARC+QCR turbulence 
model predicts wing stall after C; may at « = 12°, and there is a break in the C, curve at a = 13° from the loss 
of lift, due to inboard flow separation. The break in the C, curve occurs at higher angle of attack (a = 16°) 
with the two-equation BSL turbulence model, and the loss of lift is more severe from the effect of both 
inboard and outboard flow separation. Figure 10 shows skin friction coefficient for a = 13° (SARC+QCR 
turbulence model) and a = 16° (BSL turbulence model) in the x direction. A cut-off limiting value is 
specified to remove reversed, separated flow from the view of the wing surface (shown in gray). Both the 
SARC+QCR and BSL turbulence models predict wide flow separation on the inboard, upper wing surface 
that extends from the leading edge to the trailing edge. Additionally, the BSL turbulence model has 
extensive outboard flow separation between the flap and the tip nacelle (aileron region) (fig. 10(b)). Figure 
11 shows the span locations (fig. 11(a)) and C,, cuts of y = 150 in., 160 in., 170 in., and 180 in. on the wing 
comparing the SARC+QCR (a= 13°) and BSL (a@ = 16°) turbulence models at the break in the C, curve. 
The pressure coefficient distributions are similar between two solutions in the blown-wing region of y = 150 
in. and y = 160 in., but the BSL solution has a higher suction peak on the upper surface leading edge because 
of the higher angle of attack (figs. 11(b) and 11(c)). The flow separation predicted on the upper surface of 
the wing is outboard of the blown-wing flow for both turbulence models. For the SARC+QCR turbulence 
model at a = 13°, the separation upstream of the aileron is almost at a constant chord location between 160 
in. < y < 180 in. At station y = 170 in., the SARC+QCR separation line is at approximately 68 percent chord 
(fig. 11(d)) and at station y = 180 in., the flow is again attached to the trailing edge. For the BSL turbulence 
model at a = 16°, the separation region extends past station y = 180 in. (fig. 11(e)) and out to the tip nacelle 
(fig. 10(b)). The separation location varies with span from 77 percent chord at y = 160 in. (fig. 11(c)), to 59 
percent chord at y = 170 in. (fig. 11(d)), and to 13 percent chord at y = 180 in. (fig. 12(e)). In summary, the 
maximum lift coefficient is higher, at higher angles of attack for BSL compared to SARC+QCR (fig. 9), but 
the loss of lift after stall is greater with BSL (AC, = 0.34) due to the addition of massive outboard flow 
separation added to the inboard flow separation present with SARC+QCR (AC, = 0.17). 

As mentioned previously, CFD grids were made for both the cruise and the high-lift wing (30° flap 
configuration), with and without the HL nacelles (fig. 3) to determine the impact of the high-lift nacelle 
geometry on the wing lift and drag. The impacts of the high-lift nacelles on the C, of the unblown wing are 
shown in figure 12 for FUN3D with the SARC+QCR turbulence model. In general, the nacelles reduce the 
lift across the range of angle of attack and the lift curve slopes are about the same for both the cruise and the 
high-lift wing. The unpowered cruise wing with high-lift nacelles has a AC, decrement of 0.10 to 0.15 over 
the computed angle-of-attack range, with the angle-of-attack at stall being the same for both cases (fig. 
12(a)). Similarly, the unpowered high-lift wing with HLN configuration has a AC; decrement of 0.05 to 
0.11 (fig. 12(b)) compared to the high-lift wing without HLN. The C, decrement due to HLN diminishes as 
C,,max 18 approached. Figure 12 also confirms the benefit of the flap when comparing the lift coefficients 
between the unblown cruise and high-lift wings. Figure 13 shows a comparison of the spanloading between 
the unblown cruise wing with and without HLN for 150 KTAS @ = 10°. The smooth blue curve 
(s4.1g1m233a0p0.sect_forces) is the c; distribution for the cruise wing without HLN and the red curve 
(s4.1g2m233a0p0.sect_forces) is the c; distribution for the cruise wing with HLN. The effects of the nacelles 
are clear with large slope changes at each nacelle location. The lift decrement is AC, = 0.09 due to the 
presence of the HLN, and is obvious in figure 13 by the reduction in c; from the blue to the red curve, with 
the largest reduction occurring between 20 < y < 180. Figure 14 shows a comparison of the spanloading 
between the unblown high-lift wing with and without HLN for a = 10°. The effect of HLN on the spanwise 
c; of the unblown high-lift wing is slightly smaller (A C, = 0.08) than the effect on spanwise c; for the cruise 
wing, and the reduction in c; is limited to 25 < y < 170. 

The impacts of the high-lift nacelles on the C, of the blown, high-lift wing are shown in figure 15 for 
FUN3D with the BSL turbulence model. The nacelles reduce the lift curve slope and delays the stall angle of 
attack for the blown high-lift wing. The powered high-lift wing with high-lift nacelles (fig. 3(b)) has a 
decrement to lift of AC; = -0.13 at a = 8° to AC, = -0.25 at a = 12° compared to the powered high-lift wing 
without HLN (fig. 3(d)). Table 4 shows the component contributions to lift coefficient for the blown high- 
lift wing with and without the HLN at 58 KTAS and a@ = 12°. The lift from the wing and flap for the high-lift 
wing without the HLN configuration is 0.349 higher than the configuration with HLN, but the contribution 
of the HLN (0.096) decreases the difference between the two configurations to AC; = -0.25. Figure 16 shows 
the comparison of pressure coefficient between the blown high-lift wing with and without HLN at a few 
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span stations for 58 KTAS and a = 12°. Comparisons were made from y = 0 in. to y = 180 in. by 10 inch 
increments, and the comparisons with little difference in C, are not included, in particular inboard and 
outboard of the blown-flow region (downstream of the propeller slipstreams). There is more suction on the 
wing leading edge and on the upper flap surface for the blown, high-lift wing without HLN, than with HLN, 
as seen by the pressure coefficient comparisons for various stations in figure 16. Computational efforts are in 
progress to reduce the impacts of the nacelle pylon geometry on the wing lift and drag coefficients. 

Figure 17 shows the effect of angle of attack on skin friction and lift coefficients for the blown high-lift 
wing without HLN at 58 KTAS for solutions from FUN3D with the BSL turbulence model. The blown, 
high-lift wing (30° flap configuration) without high-lift nacelles has C, max = 4.642 at a = 12° (fig. 15) and 
some flow separation inboard of the flap and outboard of the flap at the aileron (fig. 17(a)). The lift 
coefficient drops to the same value as the HLN configuration for a = 13° and a = 14° (fig. 15) when the 
inboard wing stalls (fig. 17(b)), and has a second drop in lift coefficient at a = 15° (fig. 15) when the 
outboard wing stalls (fig. 17(c)). Flow in the blown region of the wing (excluding the flap element) remains 
attached over the complete upper surface of the wing for both turbulence models over the range of angle of 
attack. Figure 18 further illustrates the changes in flow separation that occur with angle of attack by 
reviewing the impact on spanwise c; for the blown high-lift wing without HLN at 58 KTAS for FUN3D 
solutions with the BSL turbulence model. At @ = 12° (blue line), there is some inboard and outboard flow 
separation and the c; is lower in these regions (0 < y < 20 and 160 < y < 180) than in the blown region (24 < 
y < 160). At a = 13° (red line), the c; within the unblown region of 0 < y < 24 and the blown region of 24 < 
y < 46 drop substantially as the flow separation moves upstream toward the leading edge (fig. 17(b)). At a = 
15° (green line), the c; within the unblown region of 0 < y < 24 and the blown region of 24 < y < 46 drops 
even more than for a = 13°, and the c; within the unblown region of 160 < y < 180 is now reduced as the 
outboard flow separation moves upstream toward the leading edge (fig. 17(c)). 


V. Conclusions 

The X-57 Maxwell was designed with a small, high aspect ratio cruise wing yielding a high cruise lift 
coefficient (C, = 0.75) at a = 0°. The cruise propulsors at the wingtip rotate counter to the wingtip vortex 
and reduce induced drag by 7.5 percent at a = 0.6°. The unblown C;, ma, of the high-lift wing is 2.439. When 
the DEP system is operating at 164.4 hp the conservative estimate of effective Cy ma, 18 4.202 at a = 12° 
with the blown, high-lift wing. A second, more aggressive estimate of effective C, ma, is 4.444 from the 
two-equation BSL turbulence model. Therefore, the conservative estimate of the lift augmentation from the 
high-lift, distributed electric propulsion system is 1.7, calculated from the ratio of the DEP blown C; max = 
4.202 to the unblown C; max = 2.439. The computational results from two different turbulence models 
exceed the goal of Cy max = 3.95, and indicates that the stall speed goal can be met. 

The CFD solutions were used to assess the computed drag versus the estimated drag needed to achieve 
the design cruise speed. The resulting computed wing drag coefficient of Cp = 0.02191 is ACp = 0.00076 
above the estimated target wing drag coefficient of Cp = 0.02115. However, the drag overage of 7.6 drag 
counts is well under the included drag margin of 74.9 drag counts. Therefore, it is estimated that the design 
cruise goal should be achievable by only using about 10% of the originally planned drag margin. 

The effect of the high-lift nacelles act to reduce lift coefficient by a maximum of 0.15 on the unblown 
cruise and high-lift wings, and reduce lift by 0.25 on the blown high-lift wing. 
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Table 1. The X-57 Maxwell sizing study drag component breakdown and resulting required 
aircraft drag coefficient. 


Total Airplane Drag for Mod IV Force - Newtons Force - Ib 
Margin 133 29.9 
Interference 10.8 2.43 
Induced 165 37.09 
Wing Friction 65.7 14.77 
Wing Profile 28.2 6.34 
Tail Friction 33.7 7.58 
Tail Profile 5.69 1.28 
High-Lift Nacelles 83.1 18.68 
Cruise Nacelles 33.6 7.5 
Fuselage 404 90.82 
Total 962.79 216.44 


Cp =D/(q*S) = D/ (0.5 * p * V?*S) 
p at 8,000 ft 0.00187 slugs/ft® 

Vv 150 KTAS 

Vv 253.17 ft/s 

S 66.67 ft? 

q 59.87 Ib/ft 
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Table 2. The X-57 Maxwell sizing study estimate of wing drag buildup. 


Force - Newtons Force — Ib 
Induced 165 37.09 
Wing Friction 65.7 14.77 
Wing Profile 28.2 6.34 
Cruise Nacelles 33.6 7.55 
Subtotal (no HLN) 65.76 
High-Lift Nacelles 83.1 


Subtotal (with HLN) 


=D/(q*S)=D/ (0.5 p * V2* 8) 


p at 8,000 ft 0.00187 slugs/ft® 
Vv 150 KTAS 
Vv 253.17 ft/s 
Ss 66.67 ft? 
q 59.87 Ib/ft? 


Table 3. The drag coefficients for separate components from fully turbulent and transitional 
flow solutions of FUN3D and USMS3D at unpowered cruise, at Re = 2.83 million, 150 KTAS, 
a = 0.6°. 


Cp total Cpwing Corn Cp ALN Cp pylons AC wing 
FUN3D SST 0.03282 0.02538 0.00489 | 0.00270 | -0.00015 
FUN3D SARC+QCR 0.03254 0.02496 0.00506 | 0.00256 | -0.00004 
USM3D SA 0.03349 0.02637 0.00558 | 0.00177 | -0.00023 
USM3D SST 0.03557 0.02750 0.00539 | 0.00284 | -0.00017 
USM3D LM Transition | 0.02938 0.02353 0.00545 | 0.00077 | -0.00038 -0.00397 


Table 4. The comparison of component contributions to lift coefficient for the high-lift wing 


with and without HLN, FUN3D BSL, at Re = 1.33 million, 58 KTAS, a = 12°. 
Configuration Cheff Crewing Ci flap Churn CHiN Ch props 
with HLN 4.410 3.517 0.404 0.061 0.096 0.331 
without HLN 4.664 3.818 0.452 0.063 0.000 0.331 
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Figure 1. The original Tecnam P2006T aircraft. 


7 ee 


Figure 2. The X-57 Maxwell DEP aircraft. The Tecnam P2006T fuselage and tail with the 
DEP wing system that includes the wingtip propulsors and the DEP high-lift motors. 
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(a) Cruise wing with tip nacelles, HLN (blue), and black alignment line to emphasize the 
staggered pattern of the HLN. 


— | 


x 
(b) High-lift wing with tip nacelles, 30° flap (red), and HLN (blue). 


¥ 
| } 
Xx 
(c) Cruise wing with tip nacelles, no HLN. 


i " 
x 
(d) High-lift wing with tip nacelles and 30° flap (red), no HLN. 


Figure 3. The wing configurations for CFD analysis of the X-57 Maxwell airplane. 
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BB margin 133 N 

MMB interference 10.8N 
MEE induced 165 N 
[wing friction 65.7 N 
[wing profile 28.2 N 
|_Jtail friction 33.7 N 

J tail profile 5.69 N 

MB high-lift nacelles 83.1 N 
ME cruise nacelles 33.6 N 
MMMM fuselage 404 N 


2.0 
Celt = 
0.5 a @]/No Power 
a= —&— Cruise Power 
0.0 


4 2 0 2 4 6 8 10 12 14 #16 = 18 
angle of attack, deg. 


(a) Effective Lift Coefficient 


e@=])_No Power 


Ss Se <# Cruise Power 


4 6 8 10 12 14 
angle of attack, deg. 
(b) Drag Coefficient 
Figure 5. The lift and drag coefficient for the cruise wing with high-lift and wingtip nacelles, 
150 KTAS, M = 0.233, Re = 2.83 million, h = 8000 ft, and T = 30.5°F. Cruise power propeller 
conditions of 2250 RPM, 117.38 hp, and total thrust of 230 Ibf. FUN3D SARC+QCR. 
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upper surface lower surface 


a Cc 
0.0000 0.0007 0.0014 0.0021 0.0029 0.0036 0.0043 0.0050 ¥ 0.0000 0.0007 0.0014 0.0021 0.0029 0.0036 0.0043 0.0050 


(a) USM3D LM 


upper surface lower surface 


C, 
f 
* 0.0000 0.0007 0.0014 0.0021 0.0029 0.0036 0.0043 0.0050 10:0000':0.0007" 0.0084 ::0.0021" 0.0028 °0.0096::0.0043:'0.0050 


(b) USM3D SST 


upper surface lower surface 


C, 0.0000 0.0007 0.0014 0.0021 0.0029 0.0036 0.0043 0.0050 c, 0.0000 0.0007 0.0014 0.0021 0.0029 0.0036 0.0043 0.0050 


(c) USM3D SA 


upper surface lower surface 
C, c, 0.0000 0.0007 0.0014 0.0021 0.0029 0.0036 0.0043 0.0050 
f 0.0000 0.0007 0.0014 0.0021 0.0029 0.0036 0.0043 0.0050 = " 


(d) FUN3D SST 


upper surface lower surface 


C, 


0.0000 0.0007 0.0014 0.0021 0.0029 0.0036 0.0043 0.0050 f 0.0000 0.0007 0.0014 0.0021 0.0029 0.0036 0.0043 0.0050 


(e) FUN3D SARC+QCR 

Figure 6. The comparison of skin friction coefficient between transitional flow (LM) and fully 
turbulent flow (SST, SA, SARC+QCR) on the unpowered cruise wing with high-lift and 
wingtip nacelles at 150 KTAS, M = 0.233, Re = 2.83 million, h = 8000 ft, T = 30.5°F, and a = 0.6°. 
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Figure 7. The lift coefficient for the high-lift wing with high-lift and wingtip nacelles, 58 
KTAS, M=0.0878, Re=1.33 million, h=0 ft, and T=59°F. Blown, high-lift wing power 
conditions of 4548 RPM, 164.4 hp (13.7 hp/prop), and total thrust of 596.4 Ibf. FUN3D 
SARC+QCR. 


i i I i 
—s4.1g4m088a10p0.sect_forces 
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Figure 8. The sectional lift coefficient for the unblown (blue) and the blown (red) high-lift 
wing with high-lift and wingtip nacelles, 58 KTAS, M=0.0878, Re=1.33 million, h=0 ft, and 
T=59°F. Blown, high-lift wing power conditions of 4548 RPM, 164.4 hp (13.7 hp/prop), and 
total thrust of 596.4 Ibf. FUN3D SARC+QCR. 
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angle of attack, deg. 


Figure 9. The comparison of effective lift coefficient between FUN3D turbulence models for 
the blown, high-lift wing, with high-lift and wingtip nacelles, 58 KTAS, M=0.0878, Re=1.33 
million, h=0 ft, and T=59°F. Blown, high-lift wing power conditions of 4548 RPM, 164.4 hp 
(13.7 hp/prop), and total thrust of 596.4 Ibf. 


. Skin Friction Coefficient x direction 
x 0.004 0.008 0.012 0.016 0.02 0.024 0.028 0.032 0.036 0.04 


(a) SARC+QCR turbulence model, a = 13° 


Skin Friction Coefficient x direction 


i - 
x 0.004 0.008 0.012 0.016 0.02 0.024 0.028 0.032 0.036 0.04 


(b) BSL turbulence model, a = 16° 


Figure 10. The skin friction coefficient from FUN3D for the blown, high-lift wing with 
wingtip nacelles and HLN, 58 KTAS, M=0.0878, Re=1.33 million, h=0 ft, and T=59°F. Blown, 
high-lift wing power conditions of 4548 RPM, 164.4 hp (13.7 hp/prop), and total thrust of 
596.4 lbf. 


American Institute of Aeronautics and Astronautics 


AIAA Aviation Forum 
Denver, CO 


7150 160 170 180 


(a) Span locations y = 150 in. to y = 180 in. 


y=150 in, o=13°, SARC+QCR 
y=150 in, o=16°, BSL 


4 \ 
154 158 162 166 170 174 178 
x, inches 


(b) y = 150 in. 


y=160 in, a=13°, SARC+QCR 
y=160 in, o=16°, BSL 


4 \ 
154 158 162 166 170 174 178 
x, inches 


(c) y= 160 in. 
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y=170 in, a=13°, SARC+QCR 
y=170 in, 0=16°, BSL 
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154 158 162 166 170 174 


x, inches 


(d) y=170 in. 


y=180 in, a=13°, SARC+QCR 
y=180 in, o=16°, BSL 
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4 
154 158 162 166 170 174 
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(e) y= 180 in. 


Figure 11. The comparison of pressure coefficient between solutions at the C; break a for 
FUNS3D solutions with SARC+QCR (a@=13°) and BSL (a=16°) at outboard span stations for 
the blown, high-lift wing, wingtip nacelles and high-lift nacelles at 58 KTAS, T=59°F, h=0 ft, 
and Re=1.33 million. Blown, high-lift wing power conditions of 4548 RPM, 164.4 hp (13.7 


hp/prop), and total thrust of 596.4 Ibf. 
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@lHLN 


angle of attack, deg. 


(a) Cruise wing at 150 KTAS, Re=2.83 million, h=8000 ft, T=30.5°F 


e=¢=No HLN 


@lfHLN 


angle of attack, deg. 


(b) 30° flap configuration at 58 KTAS, Re=1.33 million, h=0 ft, T=59°F 


Figure 12. The impact of the high-lift nacelles on lift coefficient for the unblown wing, 
FUN3D SARC+QCR. 
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Figure 13. The impact of the high-lift nacelles on the sectional lift coefficient for the unblown 
cruise wing, at 150 KTAS, @ = 10°, Re=2.83 million, h=8000 ft, and T=30.5°F. FUN3D 
SARC+QCR. (Blue, without HLN; Red, with HLN) 


3.2 T T T T T I I i I 
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Figure 14. The impact of the high-lift nacelles on sectional lift coefficient for the unblown high- 
lift wing (30° flap), at 58 KTAS, a = 10°, Re=1.33 million, h=0 ft, and T=59°F. FUN3D 
SARC+QCR. (Blue, without HLN; Red, with HLN) 
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Figure 15. The impact of the high-lift nacelles on effective lift coefficient from FUN3D at 58 


KTAS, Re=1.33 million, for the blown, high-lift wing with power conditions of 4548 RPM, 
164.4 hp (13.7 hp/prop), and total thrust of 596.4 Ibf. 
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(a) y =30 in. (c) y = 100 in. (e) y = 140 in. 


160 170 760 170 1 760 170 
x, Inches x, Inches x, Inches 
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Figure 16. The comparison of pressure coefficient between the blown, high-lift wing with and 
without high-lift nacelles, FUN3D BSL, 58 KTAS, Re=1.33 million, a = 12°. 
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(a) & = 12°, Crorf = 4.642 


Skin Friction Coefficient x direction 


x 0.004 0.008 0.012 0.016 0.02 0.024 0.028 0.032 0.036 0.04 


(b) a = 13°, Cheff = 4,451 


Skin Friction Coefficient x direction 


0.004 0.008 0.012 0.016 0.02 0.024 0.028 0.032 0.036 0.04 


(c) @=15°, Chere = 4.403 
Figure 17. The effect of angle of attack on skin friction and lift coefficients from FUN3D BSL 
for the blown, high-lift wing without high-lift nacelles at 58 KTAS, Re=1.33 million and 
power conditions of 4548 RPM, 164.4 hp (13.7 hp/prop), and total thrust of 596.4 Ibf. 


Figure 18. The spanwise sectional lift coefficient from FUN3D BSL for the blown, high-lift 
wing without high-lift nacelles for a = 12° (blue), a = 13° (red), and a = 15° (green) at 58 
KTAS, Re=1.33 million and power conditions of 4548 RPM, 164.4 hp (13.7 hp/prop), and 
total thrust of 596.4 Ibf. 
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